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Abstract 

Two existing global medium-energy nucleon-nucleus phenomenological optical model 
potentials are described and compared with experiment and with each other. The first 
Q\ , of these employs a Dirac approach (second-order reduction) that is global in projectile 

^~^ [ energy and projectile isospin and applies to the target nucleus ^*^®Pb. Here the stan- 

T^ ' dard S-V (isoscalar-scalar, isoscalar-vector) model has been extended to include the 

rj- ■ corresponding isovector components by introduction of a relativistic Lane model. The 

■ I determination of the energy range, energy dependence, and isospin dependence are 

,_^ ' discussed, as are the predictions for neutron scattering observables, and also the cor- 

relations and ambiguities found in Dirac phenomenology. The second of these employs 
►^ I a relativistic equivalent to the Schrodinger equation (including relativistic kinematics) 

1/^ I that is global in projectile energy, projectile isospin, and target (Z,A). Here, particular 

cn ■ attention is given to predictions for the integrated scattering observables - neutron 

Z^ , total cross sections and proton total reaction cross sections - and their sensitivity to 

f^ I the absorptive parts of the potential. Finally, current work is described and the influ- 

t^^ ■ ence of the nuclear bound state problem (treated in relativistic mean field theory) on 

^i^ . the Dirac scattering problem is mentioned. Spherical target nuclei are treated in the 

"^ I present work and strongly-collective target nuclei (rotational and vibrational) requiring 

,_iH ' coupled-channels approaches will be treated in a future paper. 

^ : 

S I A Global Phenomenological Dirac Potential 

. ,-t ', The potential described in this section consists of a global medium-energy nucleon-nucleus 

rS I phenomenological Dirac potential for the target nucleus ^°^Pb. The potential is global in 

c3 ■ projectile energy and projectile isospin and it was determined by least-squares adjust- 

ment of calculated scattering observables (model parameters) with respect to corresponding 
measured scattering observables for both proton and neutron scattering over a wide range 
in projectile energy. 

The Dirac equation is used in the mean field approximation by which the nucleon (meson) 
fields are replaced by their expectation values. Proton-nucleus (or neutron-nucleus) scatter- 
ing is then described using isoscalar-scalar and isoscalar-vector mean fields. Here these are 
taken, respectively, as a spherically symmetric complex Lorentz scalar potential SQ{r,E, ...) 
corresponding to the (fictitious) a meson field and a spherically symmetric complex Lorentz 
vector potential Vo{r,E, ...) (time-like component of Lorentz four-vector) corresponding to 
the u meson field, together with a spherically symmetric Coulomb potential Vc- 

However, a description of nucleon-nucleus scattering requires the explicit addition of isovector- 
scalar and isovector-vector potentials (mean fields) Si{r, E, ...) and Vi(r, E, ...), respectively. 



yielding 

^ = ^o±e^i (1) 

V = Vo±eVi (2) 

e = Af-T/A={N -Z)/A. (3) 

In these equations Si and Vi correspond to 5 meson and p meson mean fields, respectively, 

\ ,Ts{proton) = -| . 



and we use the nuclear physics isospin convention: T^{neutron) = +\ ,T^{proton) — ^ 



Equations (l)-(3) are a relativistic generalization of the Lane model 

With this scalar- vector interaction including isospin the Dirac equation becomes {h = c = 1) 

[a-p + ^{m + S}]i) = [E-V -Vc]ip (4) 

where ^/^ is a four- component Dirac spinor with upper and lower components ipu and ipi , E 
is the total energy of the scattered nucleon in the cm. frame, a and (3 are four Hermitian 
operators acting on the spin variables alone (these are related to the Dirac 7 matrices), and 
ip contains a two-component isospinor which is an eigenvector of ts appearing in 5" and V . 
A second-order reduction for the upper component ipu yields 

[p" + U, + Uso{{^ ■ L) - i{r-f))}]i,u = [{E - V,f - m>f; (5) 

where the effective central potential Uc is given by 

f/c= [2E\/ + 2m5-y2 + 52 + f/ec]/2E (6) 

the Coulomb correction term U^c (numerator) is 



and the spin-orbit term f/^o is 




[/e, = -2V,V (7) 
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2E \r E^m^S-V -V^dr 



It is worth noting that 5* and V appear both linearly and quadratically in the effective cen- 
tral potential Uc leading naturally to the "wine-bottle" shapes required to describe medium- 
energy nucleon-nucleus scattering somewhat below the transition region (where the sign of 
Uc changes) ||^. Also, the Coulomb correction and spin-orbit terms both appear naturally in 
the Dirac formalism whereas they are ad hoc in the Schrodinger formalism. 

Equation (5) is solved for the extensive ^osp]^ data set by making the following assumptions 
(due to tractability and the fact that there exists much more proton data than neutron data): 
(1) the geometries of the potentials are independent of projectile species and projectile energy, 
so that all energy dependence and isospin dependence is contained in the strengths of the 
potentials, and (2) the same geometry exists for the isoscalar and isovector components of a 
given potential. With these assumptions 

U = Uo{T,e)g{r) (9) 
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where Uq is a strength, T is the projectile kinetic energy in the laboratory system, and g{r) 
is a geometric form factor taken to be a symmetrized Woods-Saxon shape (which has a 
closed-form Fourier transform) given by 

g(r) = [1 + exp (^-^^]r^ [1 + exp (-^-^]r^ (10) 

where c and a are the radius and diffuseness parameters, respectively, and c is assumed to 
be of the usual form c = tq As with r^ constant and A the target mass number. 

Energy dependence was studied by considering p -|- ^"^^Pb scattering data only which implies 
that Uo(T, e) = Uo(T) , and six forms of Uq{T) were tested : 

f/o(T) = f/o (11) 

= Uo + aT (12) 

= Uo + aln{T) (13) 

= Uoexp{-T/a) (14) 

= Uo[l + {T/af]-' (15) 

= Uoll + iT/af]--^ (16) 

The measured proton scattering observables used in studying the energy dependence consist 
of differential elastic scattering cross sections da/dQ, analyzing powers Ay{6), spin-rotation 
functions Q{0), and total reaction cross sections an, all as a function of laboratory proton 
energy Tp for the ^°^Pb target. Considering experimental data over the range 80 to 800 
MeV the minimum values of chi-square/point/energy (data set) are shown in Fig. 1 for 
four of the six energy dependencies chosen for study [the other two, Eqs. (15) and (16), 
yielded poorer results, on average, than the linear, log, or exponential energy dependencies] . 
The figure shows that no energy dependence is inadmissable and that none of the three 
energy dependencies shown is admissable over the entire energy range shown. The results of 
reducing the proton energy range are shown in Fig. 2 for these three energy dependencies. 
Clearly, a factor ~ 5 improvement in the total chi-square/point is obtained by reducing 
the proton energy range from 80-800 MeV to 95-300 MeV. Note, however, that the total 
chi-square/point is almost equivalent for the energy range of 95-500 MeV. The remainder of 
this section addresses the smallest of these three energy ranges: 95-300 MeV. 

The isospin dependence was studied by including the n + ^°^Pb scattering data (consisting 
of neutron total cross sections cr-r as a function of laboratory neutron energy T„) with the 
proton data, for the energy range 95-300 MeV. Two energy dependencies were chosen to 
study the isospin dependence. These are the logarithmic energy dependence, Eq. (13), 
which yields the best chi-square/point of all cases studied and has an energy-independent 
isospin dependence by construction (as does the linear assumption which has a slightly worse 
chi-square/point) and the exponential energy dependence, Eq. (14), which has an energy 
dependent isospin dependence by construction. Thus, the two potential strengths tested are 
of the form 

f/o(T,e) = B±eC + aln{T) (17) 

= [B±eC]exp{-T/a) (18) 
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where B, C, and a are the constants to be determined for each of the four terms of the 
complete complex scalar-vector interaction potential. The chi-square minimization led to a 
logarithmic model that gives slightly better fits to the neutron data (xtoti^) /P'^^^'^ ^^ 0-90 
vs 0.95) and the proton data (xtoi(p) /point of 12 vs 13) than the exponential model and a 
total-chi square, Xtot/point, for combined neutron and proton data, of 12.0 for the logarithmic 
model vs 12.8 for the exponential model. The best-fit potentials for these two choices of the 
energy and isospin dependence, Eqs. (17) and (18), are given in Table I. 



Table I. Best-fit Dirac global optical potentials for nucleon plus ^osp^j 
scattering in the energy interval 95 < T < 300 MeV.Q 





Logarithmic Model 


Exponential Model 


Scalar Real 


SR = -570 T 307e + 23.1/n(T) 
ro = 1.105 a = 0.692 


SR = (-491 T 362e)exp(-T/5440) 
ro = 1.102 a = 0.700 


Scalar Imag. 


SI = 237 T 71. le - 42.0/n(T) 
ro = 1.157 a = 0.512 


SI = (52.6 T 125e)exp(-r/164.2) 
ro = 1.153 a = 0.488 


Vector Real 


VR = 532 ± 235e - 37.4/n(T) 
ro = 1.109 a = 0.664 


VR = (399 ± 287e)exp(-T/1686) 
ro = 1.105 a = 0.676 


Vector Imag. 


VI = -189 ± 54.2e + 28.9/n(r) 
ro = 1.149 a = 0.633 


VI = (-54.8 ± 60.4e)exp(-T/512.2) 
ro = 1.137 a = 0.647 



Figures 3 and 4 show fits to p + ^osp]-, (j^^a at 200 MeV using the proton-only potentials and 
the neutron-plus-proton potentials (Table I) for both the logarithmic and exponential models. 
As can be seen, the fits are quite good for both models for both input data sets. In fact, on 
the basis of this 200 MeV proton data, one cannot determine the preferred model and there 
appears to be only a slight preference for the neutron-plus-proton input data over the proton- 
only input data. However, the fits to the n -|- ^°^Pb total cross section data for the identical 
two models and identical two input data sets. Fig. 5, show only qualitative agreement in 
the case of the proton-only input data whereas quite good agreement is obtained in the case 
of the neutron-plus-proton input data. In addition, the data indicate a slight preference for 
the logarithmic model over the exponential model. Furthermore, the predictive power of the 
identical two models and identical two input data sets is tested against n + ^ospj-, (differential 
elastic cross section and analyzing power data at 155 MeV [Q that were not included in the 
input data. Figure 6 shows that both potentials give similarly good predictions, but that the 
analyzing power data clearly prefers the neutron-plus-proton input data. Thus, Figs. 3-6 
lead to the conclusions that a medium-energy phenomenological nucleon-nucleus potential 
may be the best way to proceed and (somewhat weaker) that an isoscalar logarithmic energy 
dependence and an isovector energy independence may be more physical than an exponential 
energy dependence for both isoscalar and isovector components. 

Given these conclusions the logarithmic model (Table I) was used to predict unmeasured 
neutron elastic scattering angular distributions, analyzing powers, and spin-rotation func- 
tions at 100, 200, and 300 MeV [||. These are shown in Fig. 7 for 100 MeV as are the 

^Strengths are in MeV and geometry is in fm; the upper (lower) signs refer to neutrons (protons). 
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corresponding predictions for proton scattering also using the logarithmic model. The differ- 
ences between the three observables for neutron and proton scattering, at 200 and 300 MeV 
as well as 100 MeV, were studied by also performing calculations for a "gedanken" projectile 
with potential strengths appropriate to a proton, but with the charge set to zero. The study 
concluded that (a) the shift in the first minimum of the differential cross sections is due to 
the influence of the Coulomb interaction, while the enhanced magnitude of the back-angle 
neutron cross sections results from the difference in sign of the isovector strengths, (b) the 
saturation of the neutron analyzing powers (+1.0) appears to come solely from the absence 
of the Coulomb interaction, and (c) the damping of the large-angle oscillations of the neutron 
spin-rotation functions largely arises from the difference in sign of the isovector strengths, 
although the absence of the Coulomb interaction plays some role. 

Finally, the correlations and ambiguities found in Dirac phenomenology were studied [^ for 
a single case, that of p + ^°Ca at 181 MeV. Briefly, two equivalent families of potentials 
are found, only one of which predicts the correct total reaction cross section (the measured 
value was not used in determining the best-fit parameterization), and has a just slightly 
lower x^ than that of the minimum in the other family. As one might expect, relatively large 
ambiguities are found in the imaginary strengths and they are linearly correlated. Also, the 
real geometries are particularly stable and the real strengths are also correlated, but are 
much better determined than the imaginary strengths. The point to understand is that the 
observed total reaction cross section is able to distinguish the correct Dirac phenomenological 
potential family. 

A Global Phenomenological Schrodinger Potential 

The potential described in this section consists of a global medium-energy nucleon-nucleus 
phenomenological relativistic Schrodinger potential. The potential is global in projectile en- 
ergy, projectile isospin, and target (Z,A). It employs relativistic kinematics and a relativistic 
equivalent to the Schrodinger equation obtained by appropriate reduction of the Dirac equa- 
tion for a massive energetic fermion {m, k) moving in a localized central potential V{r) taken 
as the time-like component of a Lorentz four-vector. The resultant radial equation for the 
partial wave fiip) is given by (^ = c = 1) 




^ _ f/(p) L{L + 1) 



% p 



2 



fdp) = (19) 



where p = kr, Tc is the total cm. kinetic energy, L is the orbital angular momentum, and 
U{p) is the renormalized total (nuclear plus Coulomb) optical potential 

C/(rt=7V(r).7=l+^. (20) 

Equation (19) is formally identical to the radial equation for the solution of the non- 
relativistic Schrodinger equation for the analogous scattering problem. By way of example. 
Fig. 8 shows calculations of the proton total reaction cross section for p + ^'^Al using Eq. 
(19) in three different ways for the identical potential V{r) : (1) non-relativistic (classical 

5 



kinematics and 7 = 1), (2) relativistic kinematics (and 7 = 1), and (3) relativistic equivalent 
Schrodinger (relativistic kinematics and 7 > 1). Clearly, the 7 factor becomes increasingly 
important as the projectile kinetic energy increases [Eq. (20)]. Option (3) is used in the 
remainder of this section. 

The starting point for determining this potential was the phenomenological proton optical- 
model potential of Schwandt et al. based upon differential elastic scattering cross sections 
and analyzing powers for the mass range 24 < A < 208 and proton laboratory kinetic energy 
range 80 < Tp < 180 MeV. The potential employs standard Woods-Saxon form factors. 
The goals were to extend the mass range of the potential to 12 < A < 208, to extend the 
energy range of the potential to 50 < Tp < 400 MeV, and to transform the extended proton 
potential to a neutron potential for the same mass and energy ranges. Moreover, optimal 
reproduction of the measured integrated scattering observables, the proton total reaction 
cross section a^ and the neutron total cross section aT, was the main focus of the work. 

The approach used was to (a) adjust only the parameters of the proton central absorptive 
potential to optimally reproduce the measured total reaction cross sections, (b) perform 
these adjustments allowing only small changes in the calculated da/dQ and Ay{6), and (c) 
transform the extended proton potential to the corresponding neutron potential by use of 
the Lane model [^ and accounting for the Coulomb correction. [Since the proton starting 
potential does not explicitly contain a Coulomb correction term it is assumed that the 
term is implicitly present and, therefore, that it must be subtracted from the corresponding 
neutron potential. The correction is taken as 0.4 Z/As.] The work was performed by iter- 
ative computation, that is, a generalized nonlinear least-squares adjustment algorithm was 
not used,0 for three nuclei spanning a large mass range: ^^Al, ^^Fe, and ^°^Pb. The resultant 
potential gave reasonably satisfactory predictions for both proton and neutron scattering 
observables for other target A values in the same range [^. Further iterative computations 
were performed for six additional nuclei: ^^C, ^^O, ^''Ar, ^^Br, ^°^Ag, and ^^^Ba. The nine 
total extracted values of the imaginary diffuseness parameter aj, for Region II of the po- 
tential, were then fit by an expansion in powers of As as shown in Fig. 9 0. With this 
result, the current parameterization of the potential is given in Table II. An example using 
the potential is given in Fig. 10 for the integrated observables of ^^Fe and where "Modified 
potential" refers to Table II. 

Lessons from the Construction and Use of the Two Potentials 

Several conclusions (some of them tentative) come from the work summarized above. First, 
the medium-energy phenomenological optical potential is very forgiving, just like the low- 
energy phenomenological potential. In particular, several different projectile energy depen- 
dencies appear tractable (hnear, logarithmic, exponential, . . . ) provided the total energy 
range is not excessive. Also, in a Schrodinger phenomenology, relatively small adjustments 
can be made in the parameters of the absorptive potential to improve agreement with the 
integrated observables without catastrophic consequences for the differential elastic and spin- 
dependent observables. In addition, it appears possible to obtain approximately smooth 



^For this reason the results have not been submitted for pubUcation in a refereed journal. 
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Table II. Schrodinger global optical potential for nucleon— nucleus scattering 

in the target mass range 12 < A < 208 and in the projectile 

energy range 50 < T < 400 MeV. Q 



Real Central 


Vr = 105.5[1 - 0.1625/n(T)] t 16.5[(iV - Z)/A - 
rR = 1.125 + T/10^ T<130 
rn = 1.255 , T > 130 
an = 0.675 + 3.1T/10^ 


il + Ts)iOAZ/Ai) 


Imag. Central 


Wv = 6.6 + 2.73(T - 80)/102 + 3.87(T - 80)710^ 

Wv = 7.314 + 0.0462T , T > 140 

rj = 1.65 - 2.4T/10=^ , T < 140 

r/ = 1.17, T>140 

ai = 0.27 + 2.5T/10=^ , T < 140 

aj = 0.3537 + 0.08451A^ - 0.001835A§ , T > 140 


, T < 140 


Real Spin-Orbit 


Vso = 19.0 1 - 0.166/n(T) ± 3.75[(A^ - Z)/A 

ryso = 0.920 + 0.0305A3 

ayso = 0.768 - 0.0012T , T < 140 

avso = 0.60 , T > 140 




Imag. Spin-Orbit 


Wso = 7.5[1 - 0.248/n(T); 
rvKso = 0.877 + 0.0360/13 
awso = 0.62 





energy dependencies of the various observables with piecewise (continuous and discontinuous) 
parameterizations. [Note that this is not always the case for the transmission coefficients.] 
The Schrodinger phenomenology of Table II contains such parameterizations because the 
starting potential [0 was constructed in this way. In general, however, piecewise parameter- 
ization should clearly be avoided. Second, a Dirac phenomenology may provide physically 
realistic potentials over a wider projectile energy range than a Schrodinger phenomenol- 
ogy because the effective central potential in a second-order reduction involves squares and 
cross terms of the form factors appearing. This allows, for example, a "wine-bottle" shape. 
Third, in a Schrodinger phenomenology employing Woods-Saxon form factors, it appears 
that small adjustments in the imaginary diffuseness parameter a/ can fine tune the inte- 
grated observables with a minimal impact on the other observables. Perhaps the same is 
true in Dirac phenomenology? Fourth (and last), given the extreme sparseness of experimen- 
tal medium-energy neutron scattering differential elastic and spin-dependent observables, a 
Dirac or relativistic Schrodinger phenomenological approach that is global in (1) projectile 
energy, (2) projectile isospin, and (3) target (Z,A), and uses the existing medium-energy 
proton and neutron (total cross sections) databases, appears to be a tractable way to cal- 
culate physically realistic neutron elastic scattering observables over wide ranges in energy 
and target. Whether a Dirac or relativistic Schrodinger formalism should be used has yet to 
be determined. 



■^Strengths are in MeV and geometry is in fm; the upper (lower) signs refer to neutrons (protons); ra is 
defined just below Eq. (3). 

7 



Current Work 

Because a satisfactory global medium-energy nucleon-nucleus optical potential does not yet 
exist we are continuing our work on this goal. Currently, we are addressing an energy 
range of (perhaps) 20 MeV to (perhaps) 2000 MeV and a (spherical) target mass range of 
16 to 209. The experimental database (currently over 20000 points) consists of the sets 
{da /dVt, aR, Ay, Q} for protons and {ax, some da/dVt and Ay) for neutrons. Our approach 
is to consider both relativistic Schrodinger and Dirac phenomenology with the identical 
database in a nonlinear least-squares adjustment algorithm. Piecewise parameterizations 
will be inadmissable. We will also address a microscopic Dirac approach for the same ranges 
employing proton and neutron densities from recent work on the nuclear bound state prob- 
lem using a relativistic Hartree approach [|l^]. Here, only even-even target nuclei will be 
considered. This topic is particularly exciting because the relativistic Hartree approach that 
we use is easily extended to relativistic Hartree-Fock |T^ and, perhaps more importantly, 
we have discovered that our coupling constants are mostly natural (of order unity) when our 



Lagrangian is rewritten in a form that is based upon QCD scaling and chiral symmetry 111 
and whose validity demands naturalness. 
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Figure Captions 

Fig. 1 x^/point/data set for p + ^'^'^Pb scattering in the energy range 80-800 MeV for four 
energy dependencies. 

Fig. 2 x^gj/point for p + ^°^Pb scattering in three energy ranges for three energy depen- 
dencies. 

Fig. 3 Differential elastic cross sections for p + ^"^Pb scattering at 200 MeV in the log- 
arithmic and exponential models. The solid curves are obtained in calculations that 
use the Dirac global potential (Table I) which has been determined by simultaneously 
fitting neutron and proton data. The dashed curves are from calculations using a Dirac 
global potential in which only proton data have been fit. 

Fig. 4 Spin observables for p + ^osp^^ scattering at 200 MeV in the logarithmic and ex- 
ponential models. The solid and dashed curves have the same explanation as in Fig. 
3. 

Fig. 5 Total cross sections for n + ^osp]^ scattering from 95 to 250 MeV in the logarithmic 
and exponential models. The solid and dashed curves have the same explanation as in 
Fig. 3 and the dotted curve is the prediction of a geometric black disk model. 

Fig. 6 Differential elastic cross sections and analyzing powers for n + ^osp]^ scattering at 
155 MeV in the logarithmic and exponential models. The solid and dashed curves 
have the same explanation as in Fig. 3. The calculations shown are predictions as the 
experimental data were not used in determining the potentials. 

Fig. 7 Differential elastic cross sections and spin observables for nucleon-plus-^°^Pb scat- 
tering at 100 MeV as predicted by the Dirac global logarithmic potential (Table I). 
The solid curves are the predictions for neutron scattering while the dashed curves are 
the predictions for proton scattering. 

Fig. 8 Calculations of the proton total reaction cross section for the p + ^^Al reaction using 
the Schrodinger formalism, Eq. (19), in three different approaches with the identical 
optical potential. 

Fig. 9 Empirical values of and polynomial fit to the imaginary diffuseness parameter aj for 
Region II {Tp > 140 MeV) of the Schrodinger global potential of Table II. 

Fig. 10 Comparisons of measured and calculated integral scattering observables from the 
nucleon-plus-^^Fe reaction using the Schrodinger formalism. The original potential is 
that of Schwandt et al. and the modified potential is that of Table II. 
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